Production of superpositions of coherent states in traveling optical fields with 

inefficient photon detection 
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We develop an all-optical scheme to generate superpositions of macroscopically distinguishable 
coherent states in traveling optical fields. It non-deterministically distills coherent state superpo- 
sitions (CSSs) with large amplitudes out of CSSs with small amplitudes using inefficient photon 
detection. The small CSSs required to produce CSSs with larger amplitudes are extremely well 
approximated by squeezed single photons. We discuss some remarkable features of this scheme: 
it effectively purifies mixed initial states emitted from inefficient single photon sources and boosts 
negativity of Wigner functions of quantum states. 



I. INTRODUCTION 

Schrodinger's cat paradox is a famous illustration of 
the principle of superposition in quantum theory [lj. It 
poses the question of whether a classical object on the 
macroscopic level can be in a state of quantum superpo- 
sition. The component states composing such a superpo- 
sition should give macroscopically distinct measurement 
outcomes [2, 3]. A superposition of two optical coherent 
states with sufficiently large amplitudes of a tt phase dif- 
ference is considered a realization of such a macroscopic 
superposition and often called a "Schrodinger cat state" . 

Recently, such coherent state superpositions (CSSs) in 
free propagating optical fields have been found to be use- 
ful for various applications to quantum information pro- 
cessm 

td lRHilB IS E3, El El 13 . Quantum telepor- 
tation [ilia, ESQ > quantum computation |9L llOLRLlll . en- 
tanglement purification fl2| and concentration 0, error 
correction |13j. and remote entangling have been ex- 
tensively studied with CSSs. In particular, it was shown 
that quantum computation can be realized using only lin- 
ear optics andphoton counting, given pre-arranged CSSs 
as resources [l(J[ll|. In this approach, a qubit is defined 
to be a superposition of two coherent states, and all four 
Bell states can be perfectly well discriminated by pho- 
ton counting measurements and a beam splitter |5j, |l2| . 
This enables one to construct quantum gates in a rel- 
atively simple way based on the teleportation protocol 
[Tl| . The amplitudes of coherent states for qubits and 
resource CSSs should be carefully chosen for efficiency 
of quantum information processing. The CSSs of ampli- 
tude a > 2 are required as resources for efficient quantum 
computation with simple optical networks [Tl| . 

It is known to be extremely hard to generate a free 
propagating CSS using current technology. It is well 
known that the CSS can be generated from a coherent 
state by a nonlinear interaction in a Kerr medium |l4j . 
However, Kerr nonlinearity of currently available non- 
linear media is extremely small compared with the level 
required to generate a CSS and attenuation in the media 
is not negligible [l5|. 

Some alternative methods have been studied to gen- 



erate a superposition of macroscopically distinguishable 
states based upon conditional measurements 0, 0] . A 
crucial drawback of these schemes is that highly effi- 
cient photon detection is necessary. Both Song et a/.'s 
scheme EH and Dakna et al.'s one E3 require photon 
number resolving measurements which is extremely de- 
manding using current technology. Some other schemes 
|l8j require many single photon detectors instead of one 
n-photon counting detector. Even though it is known 
that many perfect single-photon detectors enable one to 
perform nearly perfect n-photon counting, such a scheme 
would suffer a similar difficulty due to detection ineffi- 
ciency of many single photon detectors. Many perfect de- 
tectors used to produce macroscopic superpositions can 
be replaced with two n-photon Fock states and two per- 
fect detectors This employs another unavailable 
factor (two n-photon Fock states) by current technol- 
ogy. A modified scheme [2(j of Ref. E3 suggested by 
Montina can be robust to detection inefficiency under 
certain conditions where success probability is extremely 
low and amplitudes of the generated CSSs are small such 
as a < 1. None of the above schemes based on con- 
ditional measurements are currently feasible to generate 
CSSs with high fidelity, the main difficulty being the un- 
avoidable inefficiency of photon detection. 

Cavity quantum electrodynamics has been studied to 
enhance nonlinear effects to generate macroscopic super- 
positions [2]]]. Some success has been reported in creat- 
ing such su perp osition states within high Q cavities in the 
microwave |22| and optical |23| domains. However, most 
of the schemes suggested for quantum informati on p ro- 
cessing with coherent states H U H II H H ES El III E3 
require free propagating CSSs. 

Recently, it was shown that free propagating optical 
CSSs with amplitude up to a — 2.5 and fidelity F > 0.99 
can be generated with squeezed single photons and simple 
all-optical operations [24( , where neither efficient photon 
detection nor \^ nonlinear interactions are required. It 
was also found to be resilient to photon production inef- 
ficiency to some extent as its first step effectively puri- 
fies initial mixed states emitted from an inefficient single 
photon source 24] . In a more general sense, these exam- 
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pies reveal that the first excited energy eigenstates can 
be converted to a superposition of macroscopically dis- 
tinguishable states by linear operations and projections. 

In this paper, we extensively analyze the scheme in 
Ref. [24| and find that its purification effects can last 
for further steps. The non-deterministic CSS amplifica- 
tion scheme is found to boost non-classicality of quan- 
tum states: even very small amount of negativity can be 
drastically increased by this process. It is also pointed 
out that single photon source is not necessary to ob- 
tain squeezed single photons if another non-deterministic 
technique, photon subtraction |25| as demonstrated in a 
recent experiment |26l |. is employed. 

This paper is organized as follows. In Sec. II, We 
briefly define and discuss the CSS as a macroscopic super- 
position state with Schrodinger's cat paradox. In Sec. Ill, 
it is shown that a CSS with a small coherent amplitude 
(a < 1.2) and high fidelity (F > 0.99) can be deter- 
ministically generated by squeezing a single photon. The 
discussion is motivated by the approach of Ref. 0] . The 
Wigner functions of squeezed single photons and CSSs 
are analytically obtained, and they are compared to vi- 
sualize the effects of squeezing on single photons. In 
Sec. IV, we fully analyze and discuss the CSS amplifica- 
tion scheme with beam splitters, auxiliary coherent fields 
and inefficient detectors. Sec. V combines the two ideas 
from Sec. Ill and Sec. IV to produce CSSs with ampli- 
tude a > 2. Weak squeezing, beam mixing with an aux- 
iliary coherent field and photon detecting with threshold 
detectors are enough to generate a CSS with amplitude 
up to a — 2.5 and high fidelity (F > 0.99) given a sin- 
gle photon source. Pufication effects for an inefficient 
single photon source are another remarkable aspect of 
our scheme, which will be discussed in Sec. VI particu- 
larly for multiple iterations of the process. We conclude 
with some final remarks in Sec. VII. A recent experi- 
ment by Wenger et al. ^2(5] is briefly addressed from the 
viewpoint of CSS generation. We emphasize that sin- 
gle photon sources are not necessary to generate CSSs of 
a > 2 employing the photon subtraction technique with 
our amplification scheme. 

II. SUPERPOSITIONS OF COHERENT STATES 
AS MACROSCOPIC SUPERPOSITIONS - CAN 
THEY BE CALLED "SCHRODINGER CATS"? 

A CSS can be defined as 

|CSSp(a)) = N v (a)(\a) + e iv \ - a}), (1) 

where N v (a) is a normalization factor, | ± a) is a coher- 
ent state of amplitude ±a, and (p is a real local phase 
factor. The amplitude a is assumed to be real for sim- 
plicity without loss of generality. In this paper we refer 
to the magnitude of a as the size of the CSS. Note that 
CSSs such as |CSS±(a)) = N±(a)(\a) ± | -a)) are called 
even and odd CSSs respectively because the even (odd) 
CSS always contains an even (odd) number of photons. 



In Schrodinger's paradox, a classical object (cat) is in 
a superposition of two macroscopically distinguishable 
states (alive and dead) . Leggett and Garg have shown the 
incompatibility between quantum mechanical prediction 
and macroscopic realism for a macroscopic superposition 
state 2j. The same kind of discussions have been made 
by Reid to show violation of Bell's inequality when lo- 
cal realism is macroscopically defined 3]. A CSS is in a 
macroscopic superposition state when its amplitude is ap- 
propriately large. It is often referred to as a "Schrodinger 
cat state" or simply "cat state" albeit there exists some 
dispute over the term. 

There are probably two conspicuous characteristics of 
CSSs which may justify the title "Schrodinger cat states" . 
Firstly, coherent states are known as the most classical 
states among pure states. The coherent states were orig- 
inally suggested by Schrodinger as a quantum analogy of 
classical particles [27]]. A classical particle can be repre- 
sented as a point in the phase space while it is prohibited 
by the uncertainty principle for a quantum state. A co- 
herent state provides the most point-like description of a 
quantum particle in the phase space among all quantum 
states. Furthermore, the coherent states do not change 
their localized shapes as they move in a harmonic os- 
cillator potential. Their Wigner functions are positive- 
definite and their P-function exist even though they are 
delta functions 28] . 

Secondly, two coherent states are macroscopically dis- 
tinguishable when they are well separated in the phase 
space. Homodyne detection can be considered a macro- 
scopic measurement as it does not resolve individual 
quanta (photon). The error probability P e of discrimi- 
nating two coherent states |a) and | — a) by a homodyne 
detection is [2|| 

P e = Erf(V2a) -I, (2) 

where Erf(x) is the error function. The error probability 
P e corresponds to the probability of a wrong discernment 
by the homodyne detection due to the overlap between 
the two coherent states. The probability P e is extremely 
small as P e < 3.2 x 10~ 5 for a > 2. In such a case, a 
CSS in Eq. (JIJ can be considered a superposition between 
two macroscopically distinguishable states of a classical 
system. 

The first characteristic explained above could be more 
or less weaker than the second one as a justification 
for "Schrodinger cat states" being an alternative title of 
CSSs. The coherent states are, of course, still far from 
typical classical objects. It was shown that quantum key 
distribution using coherent states and homodyne mea- 
surements is secure against any individual eavesdropping 
attack 30]. It has also been argued that a weak mea- 
surement of the squared quadrature observable may yield 
negative values for coherent states [3l| . 

The virtual cat in Schrodinger's paradox is, to be more 
precise, entangled with a microscopic quantum object 
while a CSS in Eq. Q is in a single mode superposition. 



3 



However, an entangled coherent state 

|ECS) oc \a)\j3) + e iv \-a)\-P) 



(3) 



can be simply generated by dividing a CSS using a beam 
splitter. Such entanglement in Eq. of macroscopically 
distinguishable states is perhaps more closely aligned 
with Schrodinger's original concept 0. 

We have shown that the error probability of discrim- 
inating between two coherent states, \a) and | — a), is 
extremely small for a > 2, which justifies the CSS as, 
at least, a macroscopic superposition state. This value 
(a > 2) is also appropriate for quantum computation 
using optical coherent states ^(j- Therefore, we are par- 
ticularly interested in generating CSSs of a > 2 in this 
paper. 



III. GENERATION OF SMALL COHERENT 
STATE SUPERPOSITIONS 

There have been some trials to generate macroscopic 
superpositions using the optical parametric amplifier and 
a single photon source Ilal32ll . In this Section, we show 
how a previous scheme |lOl to generate macroscopic su- 
perpositions can be significantly simplified so that small 
CSSs with high fidelity can be deterministically gener- 
ated simply by squeezing single photons. 



A. Simplified generation of small coherent state 
superpositions 

The key idea presented in this Section is motivated by 
the scheme described in 0|. This scheme uses a non- 
linear coupling between two optical modes which realizes 
a quantum non-demolition (QND) measurement on the 
Y quadrature of one mode. This is done by coupling the 
mode contai ning the signal to another ancillary mode as 
described in [33]. We define the quadrature operators of 
a single mode in terms of the usual creation and annihi- 
lation operators as X — a + and Y = —i(a — (r) so 



that [X, Y] = 2i and hence 



(AX 2 ) (AY 2 ) = 1. 



(4) 



So as to avoid confusion between the two quantized 
modes of the EM field in the QND appratus, the mode 
which the QND measurement is performed is called the 
signal mode and the mode it interacts with to assist with 
the measurement is called the meter mode. Also the op- 
erators associated with the observables of these modes 
are labelled with subscript 's' for signal and 'm' for me- 
ter. The device described in [33| uses the two polarization 
modes of a single spatial mode as the signal and meter 
modes. For example, the horizontal polarization might 
contain the signal and vertical polarization the meter. 
The meter mode is usually assumed to be prepared in 
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FIG. 1: A schematic of the simplification of the CSS gen- 
eration. PC represents photon counting and HD represents 
homodyne detection, (a) Conditional production using QND 
with photon counting, (b) conditional production using QND 
with homodyne detection, and (C) deterministic production 
only by squeezing a single photon. 



the vacuum state. The two polarization modes are mixed 
by a wave-plate by an angle 8. Then two mode squeez- 
ing is performed between the two polarization modes by 
a non-linear crystal. The squeezing parameter r is 
determined by the power applied to a pump beam which 
creates a squeezed vacuum in the absence of any input 
signal. Finally the polarizations are mixed by a wave- 
plate by the same angle 8. The evolution through this 
device is unitary and can be represented by a unitary 
operator U. When the squeezing parameter and wave- 
plate mixing angle are related by tanh r = sin 28 (called 
the QND condition) then the quadrature operators trans- 
form as ^(| 



X„. 



Y s 
Y 
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Y 



(5) 



(6) 



From these relations it is possible to see that the Y s op- 
erator is left unchanged through the apparatus but the 
Y m operator is mixed with the Y s operator. This allows 
information about the Y quadrature of the signal to be 
gathered from the meter mode whilst leaving the signal 
itself undisturbed. Note that to satisfy the uncertainty 
relation from equation 0] the X quadrature of the signal 
output is not identical to its input. 

In attempting to generate CSSs, the scheme in 0] 
suggests preparing the signal (and meter) in the vacuum 
state. Then apply the QND apparatus just described and 
perform a photon number measurement on the meter. 
The signal output is only accepted when a predefined 
number of photons is registered in the measurement. It 
is shown heuristically in 16] that one would expect a CSS 
to be generated when r > 1. When an odd number of 
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photons is counted in the meter mode the output is close 
to an odd CSS and when an even number of photons is 
counted the output is close to an even CSS. 

The scheme described in [l(| relies on photon count- 
ing measurements to post-select the desired output state. 
This ability to conditionally select the output induces 
the required non-linearity. However this requires efficient 
detection and photon number measurements which are 
difficult to implement. One possible resolution to this 
would be to use detection schemes which can be made to 
perform efficiently when post-selecting the output state. 
We suggest here to use homodyne detection which might 
be performed with much higher efficiency than photon 
counting measurements. 

Homodyne measurements are effectively a measure- 
ment of the X or Y quadratures depending on the phase 
of a reference signal. The eigenvalue spectrum of these 
operators is continuous which makes post selection on a 
particular value have little meaning as that one value is 
infinitesimally small in the set of all possible values. For 
example, one could measure X and select the result for 
eigenvalue zero but one could never be sure that the re- 
sult was precisely zero. To circumvent this problem one 
can use the technique shown in [3j| which accepts events 
within a given window of possible values and rejects all 
others. 

In order to produce states similar to those in 0] by 
homodyne post-selection we suggest preparing the sig- 
nal mode in a Fock state, then post select by performing 
a homodyne measurement on the meter mode after the 
QND apparatus. The condition to accept the output of 
the signal is if the measurement was in the range (—5,5) 
where 5 is some small constant. We provide a heuristic 
description in a similar fashion to that in [l6| but per- 
form a more complete analysis on a simplification that 
naturally arises when trying to write down the output 
state. 

The heuristic description proceeds by expanding the 
signal mode photon number operator in terms of the 
quadrature operators; i.e. 
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Y 2 

St 

4 



(7) 



where the subscript s represents the fact that the oper- 
ator represents the signal mode and the subscript i for 
input. The QND apparatus leaves the Y quadrature of 
the signal unchanged by equation © so 
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Y 2 

SO 

4 



(8) 



Here the subscript o represents the operator for the out- 
put mode. Now the X quadrature of the signal trans- 
forms as 

X so = X si - 2 sinhrX ml (9) 
and the Y quadrature of the meter transforms as 

Y mo = Y mi + 2 sinh rY sl (10) 



where the subscript m represents the meter mode. We 
now require a post-selective measurement on Y mo . We 
use a semi-classical approach to complete this heuristic 
description by converting operators back into classical 
variables. In terms of the semi-classical quadrature vari- 
ables, the effect of the post-selective measurement can be 
included by setting Y mo = 0. So we can write the signal 
output as 



Y M 



1 



2 sinh ? 



-V 

-L rr. 



(11) 



From this equation we can see that after post-selection 
the Y quadrature of the signal output is a scaled form of 
the Y quadrature to the meter input. Substituting this 
expression into equation [3 and rearranging for the signal 
output X quadrature one obtains 



X so = ±\Un si + 2- 



Y n . 



2 sinh? 



-2sinhrX m i. (12) 



Here we consider n S i the semi-classical form of the sig- 
nal input photon number operator which will only take 
on integer values. The meter input is a vacuum state 
and hence the X and Y quadratures contain Gaussian 
noise. The signal input contains a definite photon num- 
ber. So the first term under the square root will act like 
a constant. The second term will have some scaled Gaus- 
sian random noise and so will the term outside the square 
root. However, the multi- valued nature of the square root 
gives the two peaked behaviour required. This completes 
our heuristic description. 

In this paper we will not perform an in depth analy- 
sis of this device in full to rigorously confirm the results 
of the heuristic argument just given. However a com- 
plete analysis has been performed which confirms this 
description |34|. Here we will analyze a simplification of 
this device which has similar functionality and will prove 
fruitful towards achieving our goal of a simplified exper- 
iment. 

One can show that the unitary operator which gener- 
ates the operator transformation equations[S]and[n]of the 
QND apparatus is 



jj _ e i2sinh(r)X m Y 3 



(13) 



Now consider this operator acting on the initial state 
where the signal is in a Fock state and the meter is in 
the vacuum state; i.e. 



U\n) s \0) T 



(14) 



Inserting two instances of the identity, one expanded over 
the eigenstates of X m and the other eigenstates of Y s so 
we have 



oo />oo 



— oo •/ — oo 



dx m / dy s e 



i2sinh(r)X m Y s 



Us) \%m ) (Vs\n) (x m |0) . 

When expressions for the inner products are substituted, 
this equation becomes 
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/OO fOO 
dx m / dy s e as ^ x ^e-y"' 2 H n {y s )e- x - ,2a \y s ) \x m ) . 
-oo J —oo 



(16) 



r 



where jV n and N a are normalization factors for the nth 
photon term and a Gaussian with standard deviation a 
respectively and H n (y s ) is the nth Hermite polynomial. 
The process of projecting onto the Y m — state can be 
included by taking the inner product of this state with 

I 



the state on the meter mode alone dY (y m = 0|. This 
leaves us with a (y m — 0\x m ) inside the integral. This 
term is e l0 ' Xm — 1. Hence the post-selected output state 
is 



/OO /'OO 
dx m / dy s e i2sinh ^ x ™y°e-y^ 2 H n (y s )e- x2 ^ 2 ° \y„) . 
-OO J —oo 



(17) 



When the terms involving x m are grouped together and 
factored by completing the square and the integration 
over x m is performed one is left with 

/°° (l + 4<Tsinh 2 (r))i,2 
dy s e a H n (y s ) \y s ) . (18) 
-oo 

T his can be a w ritten a little clearer if we set n — 
\Jl + 4(jsinh 2 (r). Hence the signal output is 

dYjK, 



Fidelity 



/ dy s e 2 H n (y 3 ) \y s ) . 

J — OO 



(19) 



This state is not normalised as the process which we have 
chosen to generate it is non-deterministic. For the special 
case of n = 1 the normalised state is 



2^ 



dy s e <Ky 2 S> ffi(y s ) \y s ) 



(20) 



and as the Hermite polynomial H 1 (y s ) is linear, a n can 
be moved from under the square root to inside the inte- 
gral to give 



dy s e 2 Hi(ny s ) \y s ) 



(21) 



which we state is just a re-scaling (or squeezing) of the 
momentum wave-function of a single photon. So the case 
of a single photon input into the QND device with the 
projective measurement described above is equivalent to 
squeezing a single photon Fock state. 



B. Squeezed single photons and ideal CSSs 

We have thus been guided to comparing a single mode 
squeezed single photon state with an odd CSS. The single 
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FIG. 2: The fidelity between an odd CSS and squeezed single 
photon. The odd CSS is extremely well approximated by 
the squeezed single photon for a small coherent amplitude, 
a<1.2. 



mode squeezing operator is 
S{r) = 



,-,t2 



(22) 



where r is the squeezing parameter and a is the annihi- 
lation operator. This operator reduces quantum noise of 
a vacuum state in the phase quadrature by a factor of 
er r . When the squeezing operator is applied to a single 
photon the resultant state can be expanded in terms of 
photon number states as 



S(r)|l) 



^ (cosh r) 2 2«n! 



(23) 



The state contains only odd photon numbers and has co- 
efficients decaying exponentially as n increases in a simi- 
lar fashion to an odd CSS. The fidelity of this state with 
an odd CSS is 



F(r,a) = |(CSS_(a)|S(r)|l)| 2 

2a 2 exp[a 2 (tanhr — 1)] 
(cosh r) 3 ( 1 — exp [— 2a 2 ] ) ' 



(24) 
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If an odd CSS of size a is desired then the fidelity is 
maximised when r satisfies 



coshr 



1 



ly/9 + 4a 2 . 
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(25) 



Fig. |21 shows the maximized fidelity on the y-axis plot- 
ted against a range of possible values for a for the de- 
sired odd CSS. Some example values are: F = 0.99999 
for amplitude a = 1/2, F = 0.9998 for a = 1/^2, and 
F = 0.997 for a = 1, where the maximizing squeezing pa- 
rameters are r = 0.083, r = 0.164 and r = 0.313 respec- 
tively. These values correspond to V = 0.85, V = 0.72 
and V = 0.53, where V is the variance of the squeezed 
quadrature variable. Firstly note that for a very close to 
zero the fidelity approaches unity. When a — > 0, r — ► 
and hence the squeezing operator S(r) approaches the 
identity transformation. An odd CSS with a very close 
to zero has a significant contribution from a single photon 
and very little from higher odd photon numbers. This is 
the reason for the high fidelity as a tends to zero. The 
fidelity remains high for a near zero as one can match the 
three photon contribution to the CSS by the squeezing 
operator whilst still being able to neglect higher order 
photon number terms. Eventually as a increases, higher 
photon numbers cannot be matched and so as a tends to 
infinity, the fidelity tends to zero. 

The role of squeezing on single photons becomes clear 
by comparing Wigner functions of the squeezed single 
photons and CSSs. The Wigner function of a squeezed 
single photon can be obtained from its characteristic 
function 



Xs(v) = Tr\s(r)\l)(l\SHr)e^-^" 



exp [ - ^2 + e - 2 ^2 }] (1 _ ^2 + e -2r^ 

(26) 



The Wigner function is then 



W s {z) 



1 



Xs(v)d 2 V 



= - exp[-2(e 2r z r 2 + e - 2r zf)}(Ae 2r z 2 + 4 e - 2r zf - 1) 

7T 



(27) 

The Wigner function of the CSS is obtained by the same 
method as 



photon a good approximation of the odd CSS. If squeez- 
ing is too large, the CSS and squeezed single photon will 
become different. 




FIG. 3: (Color online) The Wigner functions of odd CSSs 
(left) and squeezed single photons (right). The amplitudes of 
CSSs are (a) 1/V2, (b) 1 and (c) 2. The de grees of squeezing 
of squeezed single photons are (a) 0.164, (b) 0.313 and (c) 
0.853, which are chosen for maximum fidelity with CSSs. It 
is evident from the figure that only small amount of squeezing 
makes a single photon a good approximate CSS. 



IV. NON-DETERMINISTIC CSS 
AMPLIFICATION PROCESS 



W?(z) 



tt(1 ±e~ 2a ) 



e -ia Zr+e iaz 



(28) 

where W+{z) (W c (z)) is the Wigner function of the even 
(odd) CSS. The Wigner functions of odd CSSs with am- 
plitudes l/y/2, 1, 2 and the Wi gner functions of corre- 
sponding squeezed single photons are plotted in Fig.|3| It 
shows that only small amount of squeezing makes a single 



In this Section, we show that an arbitrarily large CSS 
can be produced out of arbitrarily small CSSs using the 
simple experimental set-up depicted in Fig.0J Let us first 
illustrate this procedure with a simple example. Suppose 
that one has a collection of identical small odd CSSs with 
known amplitude a*. Two of the small CSSs are selected 
r )±2 cos4azi >, and are incident onto a 50:50 beam splitter, BS1, which 
acts on two coherent states la) and \8) as 



w-w* 24 i^ + >i-^ + >- < 29 » 

The two small CSSs are then transformed at BS1 as 
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|CSS_(a j; ))a|CSS_(a 4 ))b ^ \ ) f (\V2 ai ) g + | - V2a«) fl ) - (|V2a,-)/ + I - V^)/) |0) s 

oc |0) / |CSS+(\/2a)) 9 - |CSS+(\/2a))/|0) 9 (30) 



where the normalization factors are omitted on the right 
hand side. One can then say that if one could condi- 
tion on detecting |0) s , a larger CSS with amplitude \f2cti 
would be obtained at mode /. However, the nonzero 
overlap between the vacuum and the even CSS in Eq. i|5U|) 
will make it impossible to perform unambiguous mea- 
surements. The error due to this overlap is not negligible 
because the initial amplitude on is supposed to be a small 
value. Note that if the parity of the initial CSSs are dif- 
ferent, an unambiguous conditioning is possible using an 
ideal photodetector. The two small CSSs of different par- 
ity are transformed at BS1 as 



Detector A Detector B 



|CSS_(a i ))a|CSS+(a i )) b 



BSl 



|0) / |CSS_(V2a)) 8 + |CSS_(V2a))/|0) fi 



(31) 



where the normalization factor is omitted again on the 
right hand side. Since the overlap between the vac- 
uum and the odd CSS is zero, a larger odd CSS, 
|CSS_(\/2a)) f, can be conditionally produced regardless 
of the value of on by detecting no photon at mode g. 
Even in this case, however, the resulting states of condi- 
tional measurements will be highly sensitive to detection 
inefficiency for small on. An additional step therefore 
is required to unambiguously discriminate between the 
vacuum and coherent states | ± y/2on) g with inefficient 
detectors. Another 50:50 beam splitter, BS2, mixes the 
field at mode g and an auxiliary coherent state \y2ai) c 



|BSl) />g |V2 ai ) c ^ |0) / (|2o i ) t i|0) t2 + |0> t i|2o i > t2 ) 

- (|V2oi) / + | - V2ai) /) \a t ) tl \ - a t ) t2 (32) 

where |BS1) t g represents the right hand side of Eq. I|30|) 
and the normalization factor is omitted. Finally, pho- 
todetectors A and B are set to detect photons at modes 
tl and t2. The remaining state at mode / is selected 
only when both the detectors detect any photon(s) at 
the same time. In this case, it is obvious that the right 
hand side of Eq. l|32"|) is reduced to a larger CSSs. If ei- 
ther of the detectors fails to click, the resulting state is 
discarded. 

This process can be successively applied until a CSS 
of a sufficiently large amplitude is obtained. Suppose 
that an even CSS with amplitude a > 2 is required while 
the initial amplitude of small odd CSSs is on = 1. One 
may consider an experimental setup depicted in Fig. [S] 
to obtain a CSS of the required amplitude. Here we re- 
fer to the CSS amplification process depicted in Fig. 0] 
as "CA" . Firstly, four pairs of odd CSSs (i.e. eight odd 




( Ismail CSS(a) >, > 



Ismail CSS(P) > 



FIG. 4: A schematic of the non-deterministic CSS- 
amplification process. Two small CSSs at modes a and b are 
added to produce a larger CSS at mode / by a conditional 
measurement using detectors A and B. See text for details. 




( |cssia«)>; 



FIG. 5: A schematic of successive applications of the CSS 
amplification processes. CA represents the CSS amplification 
process depicted in Fig. [I] As the first step, eight CSSs of 
initianl amplitude en are fed into the four CA processes. If 
all the eight detectors in the four CA processes click, the 
resulting states are selected for the second step, and so on. In 
this example, 14 detectors and 14 beam splitters are required 
with 7 auxiliary coherent states to distill a CSS of amplitude 
2\/2ui from smaller CSSs of amplitude on. 



CSSs) with amplitude a, = 1 should be fed into four CA 
processes simultaneously as shown in Fig. [SJ If the first 
step in Fig. is successful, i.e., all the eight detectors in 
the first four CA processes click, two pairs of even CSSs 
with amplitude \plon will be generated out of the four 
pairs of smaller CSSs of amplitude on fed into the first 
four CA processes. In the second step, two CA processes 
are performed with the two pairs of even CSSs gener- 
ated from the first step. Note that the auxiliary coherent 
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states for the second step should be |2c\!j)'s. Through 
this second stage, one pair of even CSSs of amplitude 
2cti can be gained from the two pairs of even CSSs with 
amplitude \p2cii. Finally, an even CSS with amplitude 
2\/2ai(~ 2.83) can be generated by the third step which 
is only a single CA process with an appropriate auxiliary 
state. An CSS of an arbitrarily larger amplitude can 
be produced by increasing the number of the steps from 
any smaller CSSs. Of course, the success probability will 
rapidly drop down and the required resources will expo- 
nentially increase as the number of steps increases unless 
quantum optical memory is available. 

The CA process described above can be generalized for 
arbitrarily small CSSs with known amplitudes as already 
shown in Fig.^J Suppose two small CSSs, |CSS v (a)) and 
| CSS^ (/?)), with amplitudes a and (3, The reflectivity r 



and transmitivity t of BS1 are set to r — (3/ ' \J a 2 + f3 2 
and t — a / \J a 2 + /3 2 , where the action of the beam split- 
ter is represented by 

B a , b (r,t)\a) a \p) b \ta + rP) f \-ra + tp) g . (33) 

The other beam splitter BS2 is a 50:50 beam splitter 
(r = t = 1/V2) regardless of the conditions and the 
amplitude 7 of the auxiliary coherent field is determined 
as 

7 = 2a/3/y/a 2 + f3 2 . (34) 

The two beam splitters BS1 and BS2 then transform the 
two small CSSs with the auxiliary state as 



(3 a 



\A) 



K<P+<t>)\ -A)\\ — M — 



V2 V2 



J>4> 



CSS v (a)) a |CSS (/3)) b | 7 ) 
a 2 -(3 2 



A 



>|0)h/2 7 > + e *|-*L_£)|V2 T )|0) 



r 



(35) 



= I*) 



/,H,t2 



where A = \J a 2 + (3 2 . Here, the measurement operator 
Pti,t2 can be represented as 



Pn,n = (l-|0)(0|) t i®(l-|0)(0|)t2. 



(36) 



It is then obvious from Eq. I|35|l that the resulting state 
for mode / by the "click-click" event at il and t2 becomes 
\CSS v+ ,p{A)) oc \A) + e^+^l - A), whose coherent am- 
plitude A = \/ 'a 2 + [3 2 is larger than both a and (3. The 
relative phase of the resulting CSS is the sum of the rel- 
ative phases of the input CSSs. The success probability 
Pip,<t>{ot, (3) for a single iteration of the process above is 
simply calculated as 

P^(a,p) = ($\Pa,am 

(1 _ e ~^4^)2[ 1 + cqs(^ + 0)e- 2 ( a2 +^)] 
2(1 + cos</?e- 2 « 2 )(l + cos(j)e-' 2 ^) ' 

which is plotted for a number of different combinations 
in Fig. The success probabilities depend on the type 
of CSSs (odd or even) used and it approaches 1/2 as 
the amplitudes of the initial CSSs becomes large. It is 
interesting to note that the probability P^^{a, a) for two 
identical odd CSS inputs is always larger than ~ 0.214 
regardless of the value of a as shown in Fig. HJ1 This is 
due to the fact that each odd CSS contains at least one 
photon no matter how small its amplitude is. Multiple 
iterations will rapidly reduce the success probability. For 
example, if one needs to distill a CSS of a = 2 out of 
4 CSSs of a = 1 (a = l/\/2), the success probability 



a 



FIG. 6: The success probabilities of the CSS-amplifying pro- 
cess in Fig. [1] for the input fields of two identical odd CSSs 
(solid line), two identical even CSSs (dashed line), and even 
and odd CSSs (dotted line). 



will be ~ 0.027 (4 x 10~ 4 ). If a CSS of a = 2 is desired 
out of 16 odd CSSs of a = 1/2, the success probability 
will be only 2 x 10~ 13 . However, if quantum memory is 
available one can temporarily hold the output state upon 
success waiting for the remainder of the trials to give 
a successful result. This avoids the exponential scaling 
of the overall probability of success and for the a = 2 
case just considered the average number of steps is 1731. 
The inefficiency of photodetectors will also decrease the 
success probability while it does not affect the quality of 
the obtained CSSs. 

It is worth noting that not only an arbitrary large even 
CSS but also an arbitrarily large odd CSS can be ob- 
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FIG. 7: A section of the Wigner function for an initial even 
CSS of amplitude 1/2 (solid line), the resulting CSS after the 
first iteration (dashed line), and the resulting state after the 
second iteration (dotted line) . The maximum negative values 
are shown in the figure. One can see a radical increase of 
negativity of the Wigner function. 



tained out of small odd CSSs. If a larger odd CSS needs 
to be produced, a larger even CSS obtained from a col- 
lection of initial odd CSSs and a single initial odd CSS 
can be fed into a CA process so that a larger odd CSS 
can be obtained. If the even CSS of amplitude 2^2 ob- 
tained and the initial odd cat of amplitude 1 are used as 
the two input states in Fig. QJ an odd CSS of amplitude 
3 will then be obtained. 

We have pointed out that a small CSS approaches 
a single photon while a larger CSS is a superposition 
of macroscopically distinguishable states which can be 
considered a realization of Schrodinger's paradox. This 
means that in principle, our scheme distills a superpo- 
sition of macroscopically distinguishable states from mi- 
croscopic quantum states. It also increases non-classical 
features of quantum states. The negativity of Wigner 
functions is an indicator of non-classical features of a 
quantum state. Since an even CSS approaches the vac- 
uum state as its amplitude gets smaller and the Wigner 
function of the vacuum state is positive-definite, the max- 
imum negative value of the Wigner function is small for 
a small even CSS. In this regime of a small amplitude, 
the Wigner function of an even CSS looks almost like a 
Gaussian state. Fig.|7|shows how the maximum negative 
value of the Wigner function of an even CSS increases as 
our amplification process is iterated. 



V. AMPLIFYING SQUEEZED SINGLE 
PHOTONS FOR LARGER COHERENT STATE 
SUPERPOSITIONS 

In our earlier discussions, it was shown that the fi- 
delity between a squeezed single photon and an ideal 
small cat is extremely high. Therefore, it can be con- 
jectured that a larger CSS distilled from squeezed single 
photons by our scheme will also be very close to an ideal 
CSS. In what follows, we will show that this conjecture 
is right for a < 2.5 by analytical and numerical approx- 
imations. We first choose the initial coherent amplitude 



as cti — 1/V2. The fidelity of the initial CSS, which is 
the squeezed single photon, is then F = 0.99978 for the 
appropriate squeezing parameter r = 0.163725. 

The squeezed single photon can be represented in 
terms of an ideal CSS and error components as S(r)\l) oc 
\CSS-(^)) + 5W\3) + 5&\5) + 5^ 1 7) + • • •, where the 
error terms are 



= 1/4 



j(2fc+l) = 



0.162278 fe (2fc+ 1)! - k\ 



2 fe yfc! v /(2fc + l)!(e- 1) 



(37) 



It can be simply checked that 5 (5) = 0.0129669 is the 
dominant error term and 5^ exponentially decreases 
for k > 5. The state only with S (5 \ N{\CSS^(^)) + 

<r 5 )|5)), where N is the normalization factor, will give 
a fidelity F = 0.99983 for the odd CSS \CSS-(^=)). 
In other words, the state only with the dominant error 
term can be a good approximation of the squeezed single 
photon for a weak squeezing. We therefore use 

S(r)|l) « \* i )=N i (\CSS-(l/V2))+6 i \5)) (38) 

as the initial input state, where r = 0.163725, 5i — 0.0147 
and Ni is a normalization factor. The initial fidelity 
Fi between and the ideal odd CSS is made F{ — 
0.99978, which is exactly same as the case of a squeezed 
single photon. The resulting state is obtained as 

|*(i))/,ti,ta = Pti,aB^B^\\%) a \%) h \y/^) ct (39) 
= r &ti,*2[|*(i)><*(i)]], (40) 

where Trii : t2 denotes a partial trace of modes il and i2, a 
subscript (n) indicates the number of the iterative steps 
made to amplify the CSS, and B 1:1 = B(^, ^~). 

When each of the detectors detects only one photon, 
it is straightforward to calculate the resulting state 



ti<l|ta(l|*(i)) oc |CSS+(l))-l.mi(0.828|4}-0.561|i 

(41) 

where error terms smaller than 1/3 of the dominant er- 
ror term have been discarded as they exponentially de- 
cay. Considering the normalization, the fidelity of the 
state lUTJ is calculated to be 0.99974. Note that about 
60% of all the successful simultaneous clicks at detectors 
A and B correspond to this case, where the probability 
can be calculated by P™ n ™ = (n|(m|p(„) Jn)|m). About 
30% of the successful clicks correspond to the cases that 
detector A detects two photons while detector B detects 
one photon or that A detects one while B detects two 
photons. We can make a same approximation for these 
cases and the fidelity is 0.99975. On the other hand, the 
highest overlap with a CSS of a — 1 that can be obtained 
by simply squeezing a single photon is F — 0.99711, thus 
a clear improvement has been obtained. 

In order to calculate multiple iterations we need to 
use numerical techniques. We are using coherent states 
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fidelity maximizes for a particular number of iterations. 
Fig. [S] shows the maximum possible fidelity using this 
process in (a) and the number of steps in (b) against 
the desired a in the CSSs. For example, four iterations 
starting from the initial amplitude ai = 1/2 is required 
to gain the maximum fidelity F = 0.995 for a = 2. It is 
evident from Fig. 4 that high fidelity F > 0.99, can be 
obtained up to a — 2.5. The error rate for discrimina- 
tion between coherent states with a = ±2.5 via a classical 
measurement (homodyne detection) is only 3 x 10~ 7 . 



(b) 



VI. PURIFICATION EFFECTS OF THE CSS 
AMPLIFICATION PROCESS 



FIG. 8: (a) The maximum fidelity obtained in our scheme vs 
the coherent amplitude, (b) The number of iterations which 
gives the maximum fidelity vs the coherent amplitude. The 
improvements of the fidelity are remarkable when compared 
with Fig. H 



The single photons required for our scheme could be 
generated conditionally from a down-converter [33 ■ This 
is a process (like squeezing) and does not require 
photon number resolving detection. It should be noted 
that current technology does not produce pure single 
photon states; the single photon is always in a mixture 
with the vacuum as 



(42) 



of some bounded coherent amplitude and superpositions 
there of. Provided the coherent amplitudes are not small, 
the most significant contributions to these states are Fock 
states of low number. For computations here the lowest 
thirty Fock states were used. This provides a very good 
approximation for coherent states with a < 2.5. All 29 
possible "click" events are included for all detectors. 

If one wished to create a CSS with a particular a 
with n CSS-amplification steps, then initial CSSs with 
ai = a/y/T 1 are required. As the number of steps in- 
creases the required a, decreases. When generating a 
larger CSS out of the squeezed single photon states the 



where p is the inefficiency of the photon production. 
Hence the squeezed single photon state will also be a 
mixture with a squeezed vacuum as 

pS(r) |0) <0| £t(r) + (1 -p)S(r) |1) <1| S\r). (43) 

However, an interesting aspect of our scheme is that it 
may be somewhat resilient to the photon production inef- 
ficiency because its first iteration purifies the mixed CSSs 
while amplifying them. The initial input states for the 
CSS amplification process from the imperfect single pho- 
ton source are 



Pa,b,c 



{l- P f \Sx) (Si\ ® \Si) (Si | + P 2 \S ) (Sb| ® |S ) (So| 
+ p(l-p)(\S ) (Sol <8> |Sx) (Si\ + \Si) (Si | ® |S ) (So| 

I 



® (l7>(7l) £ 



(44) 



where |S ) = S(r) |0) and |Si) = S(r)|l). Here, the 
terms with p 2 and p(l — p) are undesired error terms 
where either (or both) of the single photons is missing. 
Note that the initial amplitude is required to be small to 
produce a larger CSS with high fidelity. Provided such 
a small amplitude, input states incident onto the beam 
splitters in our experimental setup contain approximately 
only two (or slightly more than two) photons. In such 
cases the probability of simultaneous clicks at detectors 



A and B in Fig. 21 will significantly decrease when any 
of the single photons is missing. In other words, the 
undesired cases will rarely be selected for the next it- 
eration of the amplification process. We have obtained 
numerical results for the initial amplitude oti — 1/2 as 
follows by the methods that we have already explained. 
If p = 0.4, the fidelity of the initial CSS, which is a mix- 
ture with a squeezed vacuum, is F = 0.60 but it will 
become F = 0.89 by the first iteration. Thus a larger 
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FIG. 9: (a) Purity and (b) fidelity of the input state (solid 
lines), the input output state after the first iteration (dot- 
dashed lines), and of the output after the second iteration 
(dashed line) as functions of the photon production ineffi- 
ciency p. The input states are mixed squeezed single photons 
intended to be approximate odd CSSs of amplitude a = 1/2. 
The final output states should be therefore approximate CSSs 
of a = 1. 



CSS of significantly high fidelity is produced. If p = 0.25 
(p = 0.05), the fidelity of the initial CSS is F = 0.750 
(F = 0.950) but becomes F = 0.941 (F = 0.990) by the 
first iteration. 

The purification by the first iteration is directly evi- 
dent by the probability argument, but what remains to 
be shown is if this effect is strong enough for purification 
to still be achieved for multiple iterations. For a double 
iteration, four input CSSs of on = 1/2 would be required 
to obtain an output CSS of a = 1. The numerical re- 
sults of the second iteration presented in Fig. [5] show 
that the improvements of fidelity obtained by the first 
step can remain for further iterations. In Fig. Eta), the 
purity of the input state (solid line) and the purity of the 
output state after the first (dot-dashed line) and second 
(dashed line) iterations have been plotted as functions of 
the photon production inefficiency p. Note that purity 
is defined here as Tr{p 2 }. Fig. EJb) shows the fidelity 
of the input state when compared with the ideal CSS of 
a = 1/2 (solid line), the fidelity of the output state after 
the first iteration compared with the CSS of a = \j\[2 
(dot-dashed line) , the fidelity of the output state (dashed 
line) compared with the ideal CSS of a — 1. When p = 
the purity and fidelity are high as the input states them- 
selves are pure and the output fidelity is expected to be 
high. The fidelity has decreased compared with the re- 
sult after the first iteration but it is still higher than the 
initial fidelity. For example, the fidelity will change from 
0.60 to 0.89 for the first iteration and finally to 0.72 for 
the second iteration when p — 0.4. For the range of 
probabilities shown here (p £ [0,0.5]) there is always an 



improvement in purity and fidelity. 

There exists an alternate way of achieving the same 
output state using the same input CSSs but using a 
different arrangement of amplification procedures. Two 
a = 1/2 odd input CSSs could be amplified to gener- 
ate one even a = y/l/2 CSS. Then this state could be 
combined with another a = 1/2 odd CSS to create a 
a = \/3/4 odd CSS. Then finally this state could be 
combined with another a = 1/2 odd CSS to generate 
the a = 1 even CSS output. One might expect that the 
presence of the a = \/3/4 odd CSS could make some dif- 
ferences in purity and fidelity by this method. However, 
the plots for the final output state which we have ob- 
tained by the same numerical technique are identical in 
nature to those of Fig. So the purifying effects of this 
procedure do not seem to depend on the way in which 
the output state is generated. 



VII. REMARKS 

We have studied a simple all-optical scheme to gener- 
ate a linear superposition of macroscopically distinguish- 
able coherent states in a propagating optical field |24| . 
In stark contrast to all previous schemes, this scheme 
requires neither nonlinearity nor efficient photon de- 
tection to generate a superposition of macroscopically 
distinguishable states. Furthermore, it exhibits some re- 
silience to photon production inefficiency because it puri- 
fies initial mixed states. We have found that these purifi- 
cation effects can last for multiple iterations. The non- 
deterministic CSS amplification scheme has been proved 
to boost non-classicality of quantum states: even very 
small amount of negativity can be drastically increased 
by this process. This scheme non-deterministically gen- 
erates CSSs with amplitude a > 2. However, it should 
be noted that a non-deterministic CSS source is useful 
enough for quantum information processing [Tol ITU . 

In the CSS amplification process, the zero amplitude 
coherent states that occur in the detection modes in 
Eq. (|32|l may be slightly different from zero because of im- 
perfect mode matching at beam splitters. This will lead 
to a small probability of accepting the wrong state. Good 
mode matching is a requirement in any linear optical net- 
work where one wishes to measure manifestly quantum 
mechanical effects. Highly efficient mode matching of a 
single photon from parametric down conversion and a 
weak coherent state from an attenuated laser beam at a 
beam splitter has been experimentally demonstrated us- 
ing optical fibers [39l | . Such techniques could be employed 
for the implementation of our scheme. 

The dark count rate of photodetectors will affect the 
fidelity of the CSSs. Currently, highly efficient detectors 
have relatively high dark count rates while less efficient 
detectors have very low dark count rates (3^. We em- 
phasize again that our scheme does not require highly 
efficient detectors because the inefficiency of the detec- 
tors does not affect the quality of CSSs even though it 
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decreases the success probability. Silicon avalanche pho- 
todiodes operating at the visible wavelength have rela- 
tively high efficiency and a small dark count rate, which 
is preferred in our proposal. 

Once free propagaing CSSs states are generated, they 
can be detected by homodyne measurements, which can 
be highly efficient in quantum optics experiments. Inter- 
ference fringes will appear as a signature of the CSSs in 
the statistics of the photocurrent at the detectors. 

Finally, we note that there is an alternative method 
to obtain a squeezed single photon even without a sin- 
gle photon source. It is perhaps not surprising that a 
squeezed single photon can be obtained by adding a pho- 
ton to a squeezed vacuum as a^S(r)\0) = sinhr5(r)|l). 
However, an interesting observation is that a squeezed 
single photon can also be obtained by subtracting a pho- 
ton from a squeezed vacuum. This can be shown by 
applying the annihilation operator to a squeezed single 
photon: 

aS(r)\0) = cosh r5(r) 1 1). (45) 
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